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Abstract 

We have studied the ground state structure of quantum strips within the local spin-density 
approximation, for a range of electronic densities between ~ 5x10^ and 2x10^ cm^^ and several 
strengths of the lateral confining potential. The results have been used to address the conductance 
G of quantum strips. At low density, when only one subband is occupied, the system is fully 
polarized and G takes a value which is close to 0.7(2e^//i), decreasing with increasing electron 
density in agreement with experiments. At higher densities the system becomes paramagnetic 
and G takes a value near (2e^//i), showing a similar decreasing behaviour with increasing electron 
density. In both cases, the physical parameter that determines the value of the conductance is the 
ratio K/Kq of the compressibility of the system over the free one. 

PACS numbers: 73.63.-b, 73.63.Nm, 73.21.Hb 
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I. INTRODUCTION 



The conductance of quantum wires, i.e., the hnear conductivity per unit length, is a 
subject of current interest since the progress in nanostructure technology has allowed the 
fabrication of quasi one- dimensional (ID) structures. An interesting phenomenon is the 
observation of quantization of conductance G in units of 2e^//i, which reflects the number 
of active channels in the transport measurement . This result is usually explained 

by considering the allowed energy subbands of a non interacting ID electron gas, where the 
factor of 2 is due to spin degeneracy. Interaction effects have been considered for some time, 
especially in the framework of the ID Tomonaga-Luttinger model, where it is predicted that 
the conductance is renormalized to G = ^{2e^jK)^ with a parameter 7 > 1 for attractive 
interactions, 7 < 1 for repulsive interactions, and 7 = 1 for non interacting electron gas 



. However, it has been argued 



3, y, y, III I that 7 should be unity, since the 



measured conductance is determined by the non interacting electrons which are injected in 
the wire. 

The conductivity measured in the experiment is not the response to the external applied 
electric field, but to that which results adding to the applied field that caused by the induced 
charge usually called induced polarization field. In other words, the measured linear 

conductivity is the screened one, which in the random phase approximation (RPA) used in 
previous calculations coincides with the free linear conductivity. Recent measurements 3| 
in ultra low-disorder quantum wires of finite length exhibit a conductance structure close 
to G = 0.7(2e^//i) evolving continuosly with increasing electron density to G = 0.5(2e^//i), 
the value expected for an ideal spin-split subband. The structure at G = 0.7(2e^//i) was 
already observed in low-disorder quantum point contacts (quantum wires of zero length) 



by several authors 
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iV-aisorder quantum 

Q,0,Q,yQ. 



In particular, it has been argued ^| that this 
structure is a manifestation of electron-electron interactions involving spin. Theoretically, 
the structure is interpreted as some form of spontaneous spin polarization of the system 



mediated through the exchange interaction |21| 



2^,12 



25L l2fil | . The appearance of that 



structure has been also interpreted as a manifestation of a Kondo effect in quantum point 



contacts 
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29|. Other interpretations of the plateau of G at 0.5(2eV/i) are based on 



a possible Wigner crystalization of the ID electron gas at low density j30|. The evolution 
of the structure with density in finite length wires is not well understood. For finite length 
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wires, Reilly et al. Q| find additional structures in higher subbands which suggest that 
many-body effects are enhanced in longer ID wires. Long quantum wires are ideal systems 
to address one- dimensional electron transport considering electron-electron interactions and 
also in the presence of impurities, which are relevant to study the Tomonaga-Luttinger liquid 

In this paper we present a self-consistent calculation, within local spin-density approxi- 
mation (LSDA), of the ground state of quantum wires of infinite length and finite width, and 
for this reason we prefer to call them quantum strips. Their extension in the third direction 
is neglected, as in most theoretical descriptions. We apply it to obtain the screened conduc- 
tivity in the framework of linear response theory, showing that in this approximation the 
key quantity that determines the conductance of the strip is the ratio of the compressibility 
of the system over the free one, K/Kq. Within LSDA, the screened response and K are 
determined by the exchange- correlation interaction. At low densities, only one subband is 
occupied, the system is magnetized and the conductance reduces to G = {e^ /h)^ K/ Kq, 
where K is the compressibility of the spin polarized state, yielding values for G which range 
from G ~ 0.7(2e^//i) to G ~ 0.5(2e^//i) depending on the electron density. At larger densi- 
ties, two degenerated subbands are occupied, the strip is paramagnetic and the conductance 



is given by G = {2e^ /h)y K/ Kq, where K is now the compressibility of the paramagnetic 
state, yielding values above -but close- to G = {2e'^/h) and slightly density depending in the 
region where the two subbands configuration is expected to be the physically realized stable 
phase of the strip j^^. 

The system is partially magnetized when three subbands are occupied, and a paramag- 
netic state is reached again at higher density when four subbands are occupied, and so on. 
Our method allows to calculate the compressibility in all these situations -although only 
the cases of one and two occupied subbands are presented- and yields values for G that 
are in qualitative agreement with the ones found in long wires We note that a scat- 



tering matrix approach has been recently employed to calculate the conductance through a 
semi-extended barrier or well in the wire Q]. In this case, G can be expressed in terms of 
the incident electron energy E in the form G = ^J2nTn{E), where Tn{E) is the current 
transmission coefficient for an electron incident in the nth subband. In our approach, G is 
calculated in an alternative way in which the conductance is the linear current response of 
the system to an external static electric field and hence is an intrinsic property of the strip. 



3 



like for example its dielectric constant. 

The plan of this work is the following. In Sec. II we describe the model and obtain the 
phase diagram of the strip, i.e., the energy and magnetization as a function of the ID electron 
density along the strip, as well as its free response. In Sec. Ill we obtain the conductance 
by calculating the screened response to an applied electric field in the limit of zero intensity 
(linear response), showing that it is mainly determined by the compressibility of the strip. 
Some concluding remarks are given in Sec. IV, and the calculation of the density response 
function of the strip within time-dependent LSDA is presented in the Appendix. 



II. MODEL 



We consider a single, infinitely long quantum strip in the y direction built on a two- 
dimensional electron gas (2DEG) by introducing a confining potential along the x direction. 
This potential is assumed to be parabolic ^^muj^x^. In the LSDA, the single electron wave 
functions are given by the solution of the Kohn-Sham (KS) equations 



^r(r)=e,.^r(r), (1) 



where i stands for the set of quantum numbers, but spin, that characterize the two- 
dimensional (2D) single-particle wave functions, and the two-dimensional electronic density 
of the strip is 

p(r) = Ei^r(r)r 

with r = (a;,y), -q^ = 1(-1) if a =t (i), and 

^ ^ dp{r) ^ ' 9m(r) ^ ' 

with p(r) = p^(r) + p^(j) and m(r) = p^(r) — p^{v). 

The exchange-correlation energy per unit surface £xc has been constructed from the results 
on the nonpolarized and fully polarized 2DEG [2^ in the same way as in Refs. HQ, 
i.e., using the two-dimensional von Barth and Hedin [s^ prescription to interpolate between 
both regimes. In Eq. we have used effective atomic units {h = e^/e = m = 1), where e 
is the dielectric constant and m is the electron effective mass. In units of the bare electron 
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mass TUe one has m = m*me. In this system of units, the length unit is the effective Bohr 
radius = aoe/m*, and the energy unit is the effective Hartree H* = Hm*/e^. For GaAs 
we have taken e = 12.4 and m* = 0.067, which yields = 97.9 A and H* = 11.9 meV. 
Translational invariance along the y-direction imposes solutions of Eq. (Q) of the form 

1 



(3) 



where n = 0, 1, 2, 



. is the subband index. Inserting Eq. Q into Eq. one gets 

1 o loo /"/",/,/ ) 



-VI + -ujIx^ + J J dx'dy' 



where we have introduced the band-head energy e^, 
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and the 2D density, which is y independent, reads 



P<x(x) = -E\/2(eF-en,.) m^)? ■ 

TT 

n 

The ID electron density pi is obtained integrating over x: 



1 



Pi 



TT 



(4) 



(5) 



(6) 



(7) 



The Fermi energy tp fixes the number of subbands that are filled in the ground state of the 
strip. For each value of pi, it is determined by solving the KS Eqs. (jH) self- consistently 
under the condition that Eq. ((Tj) is fulfilled. 

For an infinite strip like the one considered here, the Hartree potential 

p{x') 



Vh = j j dx'dy' - 



(8) 



- x'Y + - y'Y 

is obviously divergent. As in the case of the homogenous electron gas in two or three dimen- 
sions, this requires the introduction of a neutralizing positive background. One possible way 
is to assume that this background is such that the Hartree potential is strictly cancelled out, 
and only the exchange and correlation energy terms appear in the KS equations 

0- An- 
other possibility is to introduce a positive background that only cancels out the divergency 
in the Hartree potential. With this goal in mind, let us write: 



Vh = lim j j dx'dy' e 



' Jiiy-y'). 



p{x') 



^{x-x'Y + {y-y'y 



lim 2 / dx' p{x')KQ\q{x - x')] , (9) 



where Kq is the zeroth order Bessel function of second kind 39|. Since for small arguments 



Kciis) ~ —ln(s), one gets 



Vh = lirn —2 j dx'p{x') ln[qa) — 2 j dx'p{x')ln 



(10) 



where a is an arbitrary length that in the following will be chosen equal to the confinement 

4c| . The first term of Eq. (fTIH) is divergent and is cancelled 



length a = yl/cuo, as in Ref. 
by a similar term arising from the background charges. The second term is what we take as 
effective ID Hartree potential, which is logarithmic on an arbitrary length scale 



Vm(x) 



2 / dx p{x')ln 



X — X 



(11) 



4a42|, 



This potential has been already considered in Hartree calculations of quantum wires 
and may account for possible local charge nonneutrality in systems with a finite width 
like the strips considered here. In the following we carry out Kohn-Sham calculations for 
both extreme models, i.e., when the logarithmic Hartree term Eq. (jlip is active, and when 
this term is neglected due to the assumption of a complete cancellation of the Hartree 
potential Eq. 

The energy per unit length of the strip can be calculated as 



E 

T 



OTT TT * 



X — X 



^\2{eF - en.a)Y'''' + -\l2{eF - en.a)(^n.<T + I dx S^^ix 

n,cr 

dx p{x) Vxc{x) — j dxm{x)wxc{x) + j dx dx' p{x) p{x')ln 
The energy per electron E/N is obtained dividing Eq. (fT^ by pi. 



(12) 



A. Phase diagram of the quantum strip in the Kohn-Sham framework 

We have solved the KS Eqs. for strips confined by three different values of the harmonic 
potential, namely = 2,4, and 6 meV, and for pi densities up to filling 6 subbands. 
Values of uq between 2.5 and 3.5 meV have been determined for the gate voltage close to 
the threshold of the second subband for long quantum wires using a magnetic depopulation 
technique Q|. We discuss the case in which the logarithmic Hartree potential has been 
included, which we call Hartree model, and the case in which the Hartree potential is not 
taken into account in the KS Eqs., which we call the exchange-correlation model. 
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In Figs. we have plotted E/N and the magnetization ^ = {p\ — Pi) / Pi as a function 
of the one- dimensional electron density for the Hartree model, and in Figs. EEl for the 
exchange-correlation model. The numbers along the E/N curves correspond to the number 
of occupied subbands, and the vertical lines indicate the boundary between neighbor j and 
j + 1 subband regions. 

For both models, we have found that the transition from even-to-odd number of subbands 
neighbor regions, like 2 — > 3 or 4 5, is smooth, as the changes in ^ indicate. On the 
contrary, the transition from odd-to-even number of subbands neighbor regions, like 1 — 2 
or 3 ^ 4, is abrupt, with metastability regions -not shown in the figures- extending to 
the left and right of the crossing points. Apart from the region corresponding to the first 
subband, which is fully polarized, the magnetization reached in all the other odd subband 
regions is below the maximum value one would naively expect, i.e., 1/3 for j = 3 and 1/5 for 
j = 5. This is due to the exchange- correlation energy, which lifts the degeneracy of the t, [ 
subbands. It is also worth to note that the odd subband number regions are rather narrow, 
especially in the exchange-correlation model, getting narrower as pi increases. The relevance 
of determining the boundaries of these density regions is that, in a mean field model, they 
fix the extension of the conductance plateaus, as shown in the next subsection. 

In Figs. [7| and |H1 we show the density profiles of the 2D electronic density for the 
Hartree and exchange-correlation models, respectively, and a lateral confining potential with 
ujq = A meV. The values of the ID electronic densities are indicated, and correspond to con- 
figurations with 1, 3, and 6 occupied subbands. Also shown is the local magnetization 
^ = m{x)/p{x) for the 3 subband systems. It can be seen that in this later case, the strip 
presents an edge magnetization that increases as p decreases. 

Some LSDA calculations have found that, starting from paramagnetic two subband con- 
figurations and decreseasing further the electronic density, for certain values of the wire 
width the system mav undergo a phase transition and the ground state becomes a spin- 
density wave (SDW) Using as a guide the diagram shown in Fig. 1 of Ref. 0], we 
have checked that, at low densities, for the strips considered here, their LSDA ground state 
is the spin po 
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arized one. The same happens for the narrow wires considered e.g. in Refs. 
24 . 126^ . We will come back to this point in the next Sec. when we discuss the 
conductance. It is worthwhile to note that the presence of a SDW will modify the phase 
diagram at the border between the one and two subband regions, as this band structure will 
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disappear. 



B. Kohn-Sham response of the quantum strip 

The KS calculation discussed above allows also to evaluate the KS -free, in the mean field 
sense- linear density response xoiQ,^) to a field parallel to the strip, i.e., in the y direction, 
which involves only intraband excitations: 



I r^^^.f I I Y 

^ 'i^-'-\/'i(^F-en,a) \uj + irj - kq - /2 cu + ir] + kq + qy2 J 

where 77 is a small real quantity, and hence the longitudinal conductivity associated to it 
Q: 



with the real part given by 



L q^ L 



The conductance G is obtained from Re[a{q,uj) / L] taking the uj Q and g — > limits. In 
the small q limit, for > one gets 

Im = --^8{uj - qkp ) , (16) 

where we have defined 



kY = piep - en,.) . (17) 
From Eqs. p5H16p . taking the Fourier transform, it can be easily shown that 

Thus, in the limit a; ^ and restoring Ti, we get for the conductance 

G=^El, (19) 

where the sum runs over the occupied subbands of the strip. One thus see that in the KS 
mean field model, in the density regions we have previously determined, the conductance 
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takes the values 1/2, 1, 3/2, 2, ... in units of 2e^ /h. The odd subband regions are so narrow 
that their experimental significance may be limited, especially in the more relevant case 
of full screening (exchange-correlation model). Although the KS conductance explains the 
gross experimental features, it is however unable to explain fine details, like the structures 
commented in the Introduction. In the next section we shall try to explain these findings 
incorporating interaction effects which are missing in the KS conductance. 



III. CONDUCTANCE AND SCREENED RESPONSE 



We turn now our attention to the conductivity of a quantum strip as the screened response 
to the sum of the applied external plus the local induced field. For simplicity, we address 
only the lowest density cases, namely when the quantum strip is ferromagnetic and only one 
subband is occupied, and when the strip is paramagnetic and only two subbands (spin-up 
and spin-down) are occupied, being degenerate and equally populated. These two situations 
are the more interesting ones from an experimental viewpoint. The generalization of the 
method to the case with more than two subbands is straightforward but requires a more 
elaborate calculation. 

The calculation of the conductance from the screened response is out of any ambiguity 
when the screened response corresponds to a 'fictions' neutral system with the same inter- 
action law as the screened electron-electron Coulomb interaction [l3^. This can be fully 
achieved for the exchange-correlation model previously discussed. In the case of the Hartree 
model, where local non neutrality is assumed, the calculation of the screened response is less 
clear and strongly model dependent. For this reason, in the following we develop the theory 
only for the exchange-correlation model and just make a few comments on the Hartree model 
results at the end of this Sec. 



Our starting point is the relation 



q2 



(20) 



where XsciQ-,^) is the screened response function. It is related to the linear density response 
function x(g, tu) by[3] 

X{q,^) = -, 7 r, (21) 
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where 

v{q) = 2 J dxdx'\Mx)\^ \Mx')\^Ko[q{x - x')] (22) 

is the quasi-lD Fourier transform of the Coulomb potential (see the Appendix), (j)o{x) = 
00 (x) for the ferromagnetic case, and 4>o{x) = (pl^x) = (f)o{x) for the paramagnetic case. In 
the RPA, the electron response to the total field (the sum of the external field plus the local 
induced field), Xsc{q,^), is approximated by the free response 

Xsc{q,uj) = Xo{q,^) , (23) 

yielding 

X'"-'M= ^ , . (24) 

^ -v{q)Xo{q,uj) 

To go beyond the RPA, the short-range electron-electron correlations may be taken into 
account by modifying Eq. (j23|l as 



Xsciq, ^) = , (25) 

1 + -f^g{q,uj)xo{q,uj) 



where the dynamic local field correction Q{q,u) has been introduced j45|, |46[. Eq. (^3)) is 
the most general way to express Xsc{q,^^) in terms of the free response function xo{q,^), 
and yields for x(g, tu) 



xiQ, ^) = : ^(^)7^ — - — — : • (26) 



Xo(9,^) 
l-Y[l-Giq,uj)]xoiq,uj) 

In the following we are only interested in the u = limit, so the frequency dependence of 
the local field correction is suppressed 



Xsciq, ^) = • (27) 

1 + ^G{q) xo{q,^) 



L 

In this form, the frequency dependence of ]c.^r(q, cu) comes only from Xo{q, An important 
property of the local field correction is 



lunviq) G{q) = viO) = ^ f 1 - ^) , (28) 



where K is the compressibility of the system and Kq its free value. In the situations we 
are considering, we have Kq = l/(7r^pf) for the ferromagnetic case, and Kq = 4/(7r^pf) 
for the paramagnetic case, and the compressibility K can be calculated from the standard 
-thermodynamical- expression 

1 5 / d^E/N dE/N\ , , 
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For the locally neutral system the energy functional contains only the kinetic and exchange- 
correlation pieces plus the confining potential in the x-direction. One gets 

Pi 



E/N = —pi + — / dx£^c{x) + const. , (30) 
C Pt J 



where c = 6 for the ferromagnetic, and c = 24 for the paramagnetic case. Hence 
where we have defined 

dx^\Mx)\'. (32) 
From Eqs. and (j2Z|), we write in the small q limit 



where kp = y2(ep — eo) with eo = ej and c' = 1 for the ferromagnetic case, and with 
^0 = Co = Co c' = 2 for the paramagnetic case. Hence 

1 q'^ 1 1 

^X.c(g,-) = 7c'^^^;j3(^^ = 7^Xo(g,-,4/7) , (34) 

where 

It is now possible to calculate the real part of the conductivity and thus the conduc- 
tance along the line of Eqs. ()13H19|) . using Xsc{,q,^) / L = •y xoil,^, kp/j)/ L at the place 
of xoil,^, kp)/L. One easily gets 

for the ferromagnetic case, and 

for the paramagnetic case. In both situations the ratio K/ Kq is calculated from the expres- 
sion 

K I 

^ = > (38) 



2c' T ' 



where / is defined in Eq. (I32j) with (l)o{x) = 0o(x) for the ferromagnetic case, and 
(f)o{x) = (t)l{x) = 00 (x) for the paramagnetic case. Eqs. ()35II38|) can be also obtained 
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calculating directly the screened response function in the time-dependent local spin-density 
approximation (TDLSDA), as done in the Appendix. The above derivation is however more 
general and shows that Eqs. (j36ll37|) are model independent, although the way in which we 
have calculated K is clearly model dependent. 

The ratio 7 = \J K/ Kq is plotted in Fig. IHlas a function of pi for different values of the 
frequency of the lateral confining potential ujq. One sees from this figure that when the strip is 
in the first -polarized- subband, 7 approximately ranges from 1.5 to 1 for all the confinements 
considered here, yielding for this subband a conductance which goes from G ~ 0.7(2e^//i) 
to G 0.5(2eV/i) with increasing density, in agreement with the experimental data of 
Reilly et al. |lj]. The discontinuity of 7 and K when passing from one to two subbands 
reflects the phase transition occurring in the system. After that, the paramagnetic state 
has a conductance that, starting from values a slightly larger than 2e^/h, decreases with 
increasing density to the measured value G ~ 2e^ /h. Despite that our model for an infinite 
quantum strip is obviously an oversimplification of the actual experimental device, it yields 
a qualitative agreement with measurements, especially the observed density dependence of 
the conductance. We want to comment that at values of pi lower than the ones considered 
here, 7, as calculated from Eqs. ()35H38p . becomes imaginary because the compressibility 
is negative. This result for the compressibility is well known j^l, and it appears also for 
the electron gas in two and three dimensions, reflecting the failure at low densities, of the 
jellium model in which the background of positive ions is kept indefinitely rigid. 

The phase transition between one and two bands makes room to the appearance of a 
SDW in the system, which we have not found in the calculations because of the way we 
have carried them out in practice. The instability of the paramagnetic configuration, as 
manifested by the divergency of K / Kq^ shows that the ground state of the strip must be 
different from the paramagnetic one, and cannot yet be the ferromagnetic phase, which has a 
higher energy and only becomes the ground state at lower densities. Indeed, we have checked 
that the value of the ratio C, = fic.„/ep, where is the Fermi er^ergy of the ID electrorr 
gas at the given density, introduced in Ref. [33| to characterize the one-dimensionality of 
the system, takes values that are compatible with the presence of a SDW, as shown in Fig. 
1 of that reference. To ascertain the role that a SDW might play in the conductance of a 
quantum strip is beyond the present model, which assumes the existence of a band structure 
from the start. 



12 



Finally, we discuss the effect that the non local charge neutrality implicit in the Hartree 
model of Sect. II has on the value of 7. A straightforward generalization of the calculation 
developed in the Appendix for the exchange-correlation model shows that, for the logarithmic 
Hartree potential model (jllj) . one just has to change everywhere 



1-2 dxdx'\Mx)\ \Mx')nn 



X — X 



a 



(39) 



Differently from the exchange-correlation contribution (j52|) . which is attractive, the above 
Hartree contribution is repulsive and may change the value of 7 in a relevant way. In the 
case of two occupied subbands, the effect of the Hartree contribution is to make the ratio 
K/Kq smaller than one, and so the value of 7. Quantitatively, in the density range where 
the ground state of the system is paramagnetic, the value of 7 varies from 1 to 0.9, to be 
compared with the values slightly larger than 1 obtained with the assumption of charge 
neutrality. Both models yield values of G which are close to the measured value G ^ 2e^ /h. 
We want to stress again that the results obtained when one adds a term which causes that 
charge neutrality is locally violated are very model dependent. 

IV. CONCLUSIONS 

Within LSDA, we have studied the ground state structure of quantum strips for a range 
of electronic densities and several strenghts of the lateral confining potential. As the ID 
electron density increases, the system can be in different phases characterized by different 
values of the the magnetization. These phases correspond to the filling of an increasing 
number of electronic subbands. Due to the exchange-correlation interaction, we have found 
that for odd numbers of occupied subbands larger than one, the system adquires some edge 
magnetization. The width of the non-paramagnetic density regions is rather narrow, and it 
decreases as the density increases. 

We have employed two approximations for the electron-electron interaction commonly 
employed to address the structure of quantum wires. The first one is a logarithmic Hartree 
potential successfully used in the past to address the infrared response of quantum wires in 



the RPA 







42j |. The second one supposes that the direct electron-electron interaction is 
fully screened by a jellium background, and the system is neutral at a local scale. 

We have used this second approximation to address the conductivity of quantum strips 
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when one or two subbands are occupied in the ground state, going beyond the mean field 
or random-phase approximations. We have found, in a model independent way, that the 
physical parameter that determines the value of conductance is the ratio K/ Kq of the com- 
pressibility of the system over the free one. This result has been used to obtain the con- 
ductance G of the system in the LSDA. We have found that when only one subband is 
occupied, G takes a value close to 0.7(2e^//i) and decreases as the electron density increases, 
in agreement with experiments. When two subbands are occupied and the system becomes 
paramagnetic, G takes a value near (2e^//i), showing a similar decreasing behaviour with 
increasing electron density as for one subband. 

APPENDIX 

In this Appendix we calculate the TDLSDA density response function of a quantum 
strip for the two cases -ferromagnetic and paramagnetic- discussed in Sec. Ill, and derive 
the screened response in the same approximation. We start from the time-dependent KS 
equations j3l in an external, time oscillating field A (0'''e""^* + Oe*"^*) with O = J2iLi e~^'^^\ 
in the ^/-direction along the wire: 



where pj is the jellium density. For the ferromagnetic case (p{x, y) is the single particle wave 
function of the electrons in the lowest fully polarized subband and p = p\ In the param- 
agnetic case in which electrons fill the two lowest degenerate spin-up and down subbands, 
ip{x,y) = ip^{x,y) = {p^{x,y) and p = p^ + pK Since the time oscillating field is in the 
y-direction, electron density oscillations are induced by the field only along this direction. 
We write 



where po is the density of the unperturbed ground state which we suppose to be equal to 
the jellium density for the neutral system we are considering in this Appendix, and 



{x,y,t) 




(A.l) 




(A.2) 




(A.3) 
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with 

5p {y, t) = 5p (e'^"?^-"*) + e-'^^y-'^'^) , (A.4) 

6p being a constant to be determined. The density fluctuation Eqs. flA.2|) - ()A.4|) induce 
fluctuations in the density operator F = Y.f=i e*''^* given by 



6FiO,iu) = mmm) - {0\F\0) = J dxdye''^y[pix,y,t)-po]o = Le'^'dp , (A.5) 
where L is the length of the strip. Moreover, the dynamic polarizability is given by 

x{qM = ^- (A.6) 

To determine 5p we now use that the wave function (p{x,y,t) can be factorized in a part 
depending only on x, (f)o{x), which is not affected by the external fleld, and a part which 
depends on y and time due to the external fleld, (p{y,t). Multiplying Eq. ()A.1|) on the left 
by 4>q{x) and integrating over x we get 

lv72 , , /■^.^./^./i^ p{x',y',t)-pj 



i—ip{y,t) = l--Vl + const. + J dx dx' dy' \(f)o{x 



\l{x- x'f + (y - y'f 



We then insert p{x,y,t) of Eqs. ()A.2|) - ()A.4|) into ()A.7|) and linearize the equations. This 
means writing the self-consistent KS mean fleld entering Eq. ()A.7|) . i.e. Vks[x, y, p{x, y, t)] = 
VH[x,y,pix,y,t)] +v^c[x,y,pix,y,t)], as 

OVks 



VKs[x,y,p{x,y,t)] = VKs{x,y,po) + 



6p{x,y,t) . (A.8) 

P=PO 



dp{x,y,t) 

Therefore, from Eqs. ()A.7|) and ()A.8|) we obtain: 

i^^^{y,t) = l^-]^Vl + const. + [{v{q) + I)5p + \] (e^^^^^^-*) + e'^^'^^--*)) } (y, t) , (A.9) 
where v{q) is the Fourier transform of the Hartree potential 

f (g) = 2 I dxdx\<j)Q{x)\^\<i)Q{x)\^KQ[q{x - x')] , (A. 10) 



and / is given by 



dx^\U^)\'. (A.ll) 
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This quantity represents the exchange- correlation contribution to the residual interaction. 

Deen used to describe collective 



Neglecting this term in Eq. ()A.9|) yields the RPA which has 



excitations of quantum wires by several authors 
Eq. ()A.9|) can be rewritten as 

d , ^ ( I 



3 



47 



49|. 



t-ip (y, t) = |--VJ + const. + X' [e'^'^y--'^ + e-^(<?^--*)]| (y, t) (A. 12) 



with 



X\q) = \ + [v{q) + I]5p. (A. 13) 

Equation ()A.12|) coincides with that of a non-interacting system coupled to an external time 
oscillating field, with a coupling constant A' given by Eq. ()A.13|) . For such a system, the 
density response function is the -single particle- free response xo(q'5^) we have studied in 
Sec. II. From Eq. ()A.6|) and from the analogous relation for the free response function 

Xo(g,^) = ^, (A.14) 

we obtain 

\x{q,uj) = >:{q)xo{q,uj) = L5p . (A.15) 
The solution of these equations is the TDLSDA response function 



X^^^^^(g,^) = . (A. 16) 



The RPA response function is obtained just neglecting the interaction term / in the above 
equations: 

1 - -¥xo{q,^) 



In the RPA A' = X + v{q) 6p, and electrons respond like non- interacting particles to the com- 
bined effect of the external field A (0"''e~*'^* + Oe*'^*) plus the Coulomb local polarization field 
v{q) 6p induced by the density fluctuation. Since the screened response is by definition the 
linear response of the system to the sum of the external field plus the Coulomb polarization 
field 

A x(g, uj) = [\ + 5pviq)] Xsciq, uj)=L5p (A. 18) 

( \ Xsciq, uj) /A 1n^ 
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one sees immediately that in the RPA 

Xr^(g,^)=Xo(g,^), (A.20) 

whereas in the TDLDA one gets 

Xir^'^^u;)^ ^fr"^ (A.21) 

1 - iXo(g,w) 

The response functions described above refer to electrons in the lowest subband, xily ~ 
X^{q,i^), or filling the two lowest spin-up and -down degenerate subbands, x{Qt^) — 
X^{q,^) + xKQt^)- The formalism can be easily generalized to the case in which more 
than two subbands are filled. 
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FIG. 1: Energy per electron (meV, left scale) and magnetization (right scale) in the Hartree model 
for a parabolic confinement of a;o = 2 meV as a function of the linear density (cm~^). The regions 
separated by vertical lines correspond to the indicated number of occupied subbands. For one 
single occupied subband, the system is fully polarized. 
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FIG. 2: Same as Fig. ^for ujq = 4 meV. 
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FIG. 3: Same as Fig. ^ior ujq = Q meV. 
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FIG. 4: Energy per electron (meV, bottom panel) and magnetization (top panel) in the exchange- 
correlation model for a parabolic confinement of cjq = 2 meV as a function of the linear density 
(cm~^). The regions separated by vertical lines correspond to the indicated number of occupied 
subbands. For one single occupied subband, the system is fully polarized. 
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FIG. 5: Same as Fig. 0]for ujq = 4 meV. 
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FIG. 7: Several density profiles of the 2D electronic density as a function of x (nm) for the Hartree 
model and a harmonic frequency = 4 meV (left scale) . Also shown is the local magnetization ^ 
for the 3 subband case (right scale). The values of the ID electronic densities, which correspond 
to configurations with 1, 3, and 6 occupied subbands, are indicated. 
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FIG. 8: Same as Fig. [7|for the exchange-correlation model. 
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FIG. 9: Ratio 7 = \/K/Kq as a function of the linear density pi (cm~^) for the exchange- 
correlation model and three values of the harmonic frequency uq. The lines have been drawn to 
guide the eye. 
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